Ultrasound attenuation in gap-anisotropic systems 



J. Moreno and P. Coleman 
Department of Physics and Astronomy, Rutgers University, Piscataway, NJ 08855 

(February 1, 2008) 

Transverse ultrasound attenuation provides a weakly-coupled probe of momentum current cor- 
relations in electronic systems. We develop a simple theory for the interpretation of transverse 
ultrasound attenuation coefficients in systems with nodal gap anisotropy. Applying this theory 
we show how ultrasound can delineate between extended-s and d-wave scenarios for the cuprate 
superconductors. 
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The development of gaps with nodal anisotropy is a 
recurrent property of highly correlated electron systems. 
Bulk probes, such as the NMR relaxation rate,""^ spe- 
cific heat^ and penetration depth'^''^ indicate that gap 
nodes may be present in a wide variety of strongly cor- 
related systems including heavy fermions, strong cou- 
pling and cuprate superconductors^"'' and the narrow 
gap kondo insulators CeNiSn and CeRhSb.^ However, 
with a few noted exceptions^' we have no direct infor- 
mation about the symmetry of the gap in these strongly 
correlated systems. 

A versatile, but much under-utilized tool for probing 
electronic gap nodes is the use of transverse ultrasound 
attenuation. This method has been successfully used to 
locate the gap lines and point nodes in superconduct- 
ing UPt^}^ Surprisingly, very little work has been done 
to enable the model-independent interpretation of trans- 
verse ultrasound measurements. In this paper we revisit 
this old problem, highlighting those aspects of ultrasound 
attenuation that are model-independent and relevant to 
future experiments. 

Ultrasound attenuation probes the relaxation of elec- 
tronic momentum in a model-independent fashion. This 
information is encoded in the "viscosity tensor" , a high 
symmetry tensor with very few independent components. 
Here we develop a simple theory which links these com- 
ponents to the location of the gap nodes. We illustrate 
this theory in a vein of current interest, cuprate super- 
conductivity, showing how ultrasound measurements can 
provide an unambiguous fingerprint of gap zeroes lying 
on the diagonal of the Brillouin zone.^^ 

In a metal, the phonon strain field Uij{x) couples lin- 
early to the electron stress tensor cr'-' (x) 
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This coupling is model- independent. The stress tensor 
cr*-' {x) describes the flow of electronic momentum: its di- 
vergence governs the rate of change of electronic momen- 
tum density Vjcr*-'(x) — —Pi{x). When a sound wave 
propagates through a crystal, the dissipation rate is:^^ 
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In linear response ct*- 
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iiki where 77'^'^' is the viscos- 



ity tensor. The sound attenuation coefficient is defined as 
the ratio of the time average energy dissipation to twice 
the energy flux in the wave:^^ 



a{q, u) = rj where rj = 77*^'^' 
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p is the mass density and is the speed of sound with 
wave vector q and polarization u. From the coupling (1), 
it is straight-forward to obtain the Kubo formula for the 
viscosity tensor: 



.o-(a^-'(q)a^-'(-q)) (4) 
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with q = {q,LO + i5). 

In typical ultrasonic measurements, wavelengths A are 
hundreds of microns and substantially exceed the elec- 
tronic mean-free paths (Ze), so the attenuation is safely 
in the hydrodynamic limit, le « A. In this case the mo- 
mentum dependence of the viscosity tensor is irrelevant, 
permitting us to take the limit q — > 0.^^ 

Like the elasticity tensor, the symmetry properties of 
the viscosity tensor significantly reduce the number of 
its independent components. This tensor is symmetric, 
not only in the first and second pair of indices, but also 
under the interchange of these pairs: 



ikln 



V 



V 



ikml 



V 



(5) 



Symmetry under the crystal point group further reduces 
the number of independent components. For example, 
inversion symmetry eliminates those components with 
an odd number of identical suffixes, e.g. rj^^^v = Q, 
yxyyy _ q Square or cubic environment, 90*^ rotation 
symmetry restricts the viscosity tensor to the form 

^vki ^ ^(^-*fcjj7 gugjks^ ^ + B6'^S^''6''' (6) 

with no sum on indices implied, so for a cubic or square 
lattice the ultrasound attenuation is proportional to: 
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where C = C + A and d is the dimension. We shall re- 
strict our attention to transverse ultrasound attenuation, 
for which the last term vanishes. In other words, in a cu- 
bic crystalline environment, there are only two indepen- 
dent transverse ultrasound attenuation coefficients. In a 
two dimensional hexagonal system, 60^ rotation symme- 
try means that B must also vanish. In a 3D hexagonal 
system, polynomial terms involving and Qz must be 
added to the above expression, giving 
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where we have dropped terms which vanish for transverse 
ultrasound. These strong symmetry constraints mean 
that only a few different propagation directions are re- 
quired to measure the full electronic viscosity tensor. 

We now turn to the relationship between the quasipar- 
ticle gap structure and the electron viscosity tensor. The 
electron contribution to ultrasonic attenuation only be- 
comes substantial at low temperatures in a regime where 
a quasiparticle description of the excitations is valid. For 
a fluid of quasiparticles with dispersion Ej:, the quasipar- 
ticle group velocity is = Vj^Ej^. The traceless stress 
tensor pertinent to transverse ultrasonic measurements 
is then: 
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where ibl creates a quasiparticle of momentum k and 
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is the momentum flux of a single quasiparticle. The 
quasiparticle contribution to the viscosity tensor is then 
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where / is the Fermi distribution function and is the 
relaxation time of the quasiparticle. From this relation, 
we see that a simple expression for the viscosity for trans- 
verse ultrasound is: 



) M [^fc -ujik-ql + ik- u][v^ ■ q] 
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We shall consider the situation where the temperature 
is low enough for the quasiparticles to be entirely concen- 
trated within the gap nodes of the excitation spectrum. 
To simplify our discussion, we shall assume that the gaps 
are small enough in comparison with the Fermi energy, 
that to a good approximation vp w vfPf (dEp/def), 
where and vp are respectively, the "bare" energy and 
Fermi velocity of the excitations prior to gap formation. 
Consider the case where the gap nodes are simple points 
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FIG. 1. a) The nodes are "activated", h) q or u are per- 
pendicular to the nodes, the nodes are "inactive" 



in momentum space, located at positions Po{i)- The at- 
tenuation from a given node will depend on the orienta- 
tion of the sound wavevector and polarization. If neither 
of these vectors is perpendicular to po{i) then the quasi- 
particles at the node can couple to the sound-wave. In 
this configuration, the node is "activated". (Fig. l.a) If 
however, either the wavevector direction q or the polar- 
ization u are perpendicular to the node, it is "inactive" 
and quasiparticles at the bottom of the node will not 
couple to the ultrasoimd. (Fig. l.b) In this configura- 
tion, the attenuation produced by the node is strongly 
suppressed. For most orientations of the ultrasound, the 
nodes are "activated" , and their contribution to the ul- 
trasound attenuation may be written as: 

fjA = {VFPF? Y.{u ■ Po{i)f{q ■ Po{i)fAi{T), (13) 
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is the thermal average of the inverse product of twice 
the relaxation rate Ti{E) and the quasiparticle density 
of states Ni{E) in the gap node. 

Suppose instead the gap node is "inactive" , with the 
polarization vector ii, at right-angles to Po{i), then the 
contribution to the ultrasonic attenuation contains the 
additional factor cos^(0), where cos^ = {k ■ u). In this 
case: 



fii(i) = {vFPFf{q-Po{i)fBi{T), 



(15) 
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where 



Bi{T) = cosHff{2Ni{E)Ti{E))-^ = 
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For a point-node where the size of the node grows Un- 
early with energy, Bi{T) ~ T'^Ai{T). A similar result 
holds for a line node which lies in a plane. It thus fol- 
lows that if a configuration can be found where all nodes 
are simultaneously inactive, then the ultrasonic attenua- 
tion will exhibit a relaxation rate a factor of smaller. 
Such a situation will occur for point nodes situated at 
90'^ to each-other. It will also occur for line nodes ly- 
ing in a plane, or lying in many planes that intersect at 
90°. If the nodes do not lie in such 90° configurations, 
then ultrasound attenuation will not show an anisotropic 
dependence of the powcrlaws. 

To illustrate this discussion, we now make a more de- 
tailed application to a two-dimensional model relevant to 
the cuprate superconductors. A number of recent exper- 
iments have provided strong evidence for an anisotropic 
gap with nodes at the Fermi surface"*' Superconduct- 
ing interference experiments sensitive to the phase of the 
gap function*^ support an order parameter with d-wave 
symmetry with nodes lying at 45° to the a and b axis. 
Ultrasound measurements provide a complimentary ap- 
proach which is sensitive to the precise location of the 
nodes. In particular, it offers the potential to distinguish 
between a d. 
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state and an anisotropic s-wave state 
with nodes located either side of the 45° position, and 
there are no ambiguities associated with the interpreta- 
tion of results for bilayer compounds. 

Using the symmetry arguments advanced above, if the 
wave vector q forms an angle </> with the x-axis, the trans- 
verse ultrasound attenuation a^{4i), will have the form 



- a^(0)] sin^ 2(/- + a^(0); 
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Suppose the quasiparticle excitation spectrum has the 



BCS form E^; = 
above, this leads to 

a(g, u) 
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A|. Using the results obtained 
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whore. Il{q,u) = [p ■ u){p ■ q) and ajy and tjv are the 
normal attenuation and relaxation times. 

To represent the d-wave order parameter we have cho- 
sen the function: 



Ad(6lp) = AdCos2^p 
The anisotropic s-wave is represented by: 
A,(ep) = A,$«(^p) + (A) 



(19) 
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FIG. 2. Ultrasound attenuation for a d-wave state at two 
angles: = (solid line) and = 7r/4 (dashed line). In the 
inset the ratio r = as(7r/4)/as(0) as a function of {T /Tcf . 



where: $s(6'p) = | cos(26'p)| — (2/7r). We have treated im- 
purity scattering in the standard selfconsistent T-matrix 
appr oximat ion * * . 

For a pure d-wave state (Fig. 2) the transverse ultra- 
sound attenuation is proportional to T^'^ when the angle 
is zero (the wave vector and polarization of the sound 
wave arc parallel to the symmetry axis of the crystal) and 
it is proportional to T^-^ when the angle is 7r/4. When 
there is a finite density of impurities a flat region appears 
in both attenuation coefficients (the impurities induced 
a finite density of quasiparticle states at zero temper- 
ature), but at higher temperatures the power laws are 
unchanged. Then, in accord with our earlier arguments, 
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Even considering non-resonant scattering (taking the 
cotangent of the scattering phase shift c=l) this ratio 
shows a quadratic behavior. 

The results for a pure s-wave state are very differ- 
ent (Fig. 3). We have chosen two different anisotropic 
s- waves. (A) = ;|As and (A) = 0. Unlike d-wave pair- 
ing, here the anomalous scattering off impurities is finite 
and resonant scattering does not ever develop. Conse- 
quently, the attenuation for an s-wave gap with nodes at 
45° is even more anisotropic, it is finite when the angle 
is zero and approaches zero when (j) = n/i. On the other 
hand; when the nodes of the gap function depart from 
7r/4 the ultrasound attenuation anisotropy disappears, in 
fact at temperatures close enough to T^, the attenuation 
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FIG. 3. Ultrasound attenuation for two extended s-wave 
states: As{0p) = Ao|cos20p| at angles = 
(dot-dashed line) and (p = 7r/4 (dashed line), and 
As{9p) = Ao(|cos26ip| - (2/7r)) at = (solid line) and 
= 7r/4 (dotted line). 

at = 7r/4 actually becomes bigger that at ^ = 0. In 
neither case is the ratio {n / 4:) / {0) proportional to 

ji2 

The generality of our approach lends itself naturally 

to many other strongly correlated systems with point 
nodes. One particularly interesting case is the narrow 
gap Kondo insulators, CeNiSn and CeRhSn, where 
Miyake et aP° have suggested point-nodes on the c-axis 
as an explanation of the unusual NMR relaxation rate 
and the anisotropy of the conductivity. There arc sev- 
eral other strongly correlated superconductors that de- 
serve examination, such as VsSi and UPd2Als which 
exhibit NMR relaxation rates characteristic of gap 
lines. The heavy fcrmion superconductor U Bei^ is note- 
worthy here, early ultrasound measurements^^, found no 
anisotropy, despite the clear suggestion of line nodes from 
NMR measurements.^^ This result suggests the presence 
of domains with different orientations, and it would be in- 
teresting to repeat the measurements using field-cooling 
to obtain an aligned single domain superconductor. 

This paper has emphasized the model-independent as- 
pects of ultrasound attenuation as a probe of nodal gap 
structure. This attenuation is independent of the mecha- 
nism driving superconductivity and sensitive only to the 
intrinsic symmetries of quasiparticle excitation spectrum. 
Its tensorial character permits the measurement of sev- 
eral independent components at the same time. The 
simple methods developed provide an economic way to 
extract vital information about the gap anisotropy in a 



general class of gap-anisotropic systems, and they ap- 
pear to provide a discriminating tool for elucidation of 
gap structure in the cupratcs. 
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